Excluded- Volume Effects in Tethered-Particle Experiments: Bead Size Matters 
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The tethered-particle method is a single-molecule technique that has been used to explore the 
dynamics of a variety of macromolecules of biological interest. We give a theoretical analysis of 
the particle motions in such experiments. Our analysis reveals that the proximity of the tethered 
bead to a nearby surface (the microscope slide) gives rise to a volume-exclusion effect, resulting in 
an entropic force on the molecule. This force stretches the molecule, changing its statistical prop- 
erties. In particular, the proximity of bead and surface brings about intriguing scaling relations 
between key observables (statistical moments of the bead) and parameters such as the bead size 
and contour length of the molecule. We present both approximate analytic solutions and numerical 
results for these effects in both flexible and semiflexible tethers. Finally, our results give a pre- 
cise, experimentally-testable prediction for the probability distribution of the distance between the 
polymer attachment point and the center of the mobile bead. 

PACS numbers: 82.37.Rs 82.35.Pq 36.20.Ey 87.14.Gg 
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Single-molecule biophysics has rapidly become an exper- 
imental centerpiece in the dissection of cellular machin- 
ery. This part of the biophysics repertoire often relies, in 
turn, on the use of micron-scale beads both as a reporter 
of underlying molecular motions and as the "handle" for 
grabbing these single-molecule systems. Thus, a key part 
of the theoretical infrastructure of this field is a clear un- 
derstanding of the role that these beads play in altering 
the statistical properties of the macromolecules which are 
the real target of interest in such experiments. 

Beyond interest in the in-vitro consequences of 
tethered-particle motions, many processes within the cell 
themselves involve tethering. A notable example has 
to do with vesicular trafficking, in which molecular mo- 
tors Q] carry tethered cargoes with similar Earth-like 
proportions relative to the molecular Atlases doing the 
heavy lifting. Thus, the statistical-mechanical analysis 
performed here may prove useful for understanding in- 
vivo processes, in addition to the in-vitro consequences 
that form the main motivation for the work. 

Figure^sketches the tethered particle method (TPM). 
The main idea is that a macromolecule (for example DNA 
or some protein that translocates DNA or RNA) is an- 
chored at one end to a surface, while the other end of the 
molecular complex is attached to a bead. The observed 
motion of the bead serves as a reporter of the underlying 
macromolccular motion. This technique has been used in 
a variety of settings e.g. the examination of nanometer- 
scale motions of motors like kinesin Q or RNA poly- 
merase 0, Q , protein synthesis by ribosomes [l| , exonu- 
clease translocation on DNA 0, Q rprotein mediated de- 
formation Mi and loop formation |9| in DNA, DNA hy- 
bridization and DNA motion [llL fT3 | . The main goal 
of this paper is to show how the proximity of the reporter 
bead to the surface affects the interpretation of the re- 




FIG. 1: Schematic showing the tethered particle method, (a) 
The tether is attached to a specific point on the bead; z is the 
height of this point, r denotes the position of the center of 
the bead, (b) The vector R from the attachment point to the 
bead center can rotate and is described by two angles. These 
rotations are more constrained for small values of z. Note 
that in the figure the width of DNA is not to scale, it is much 
smaller in real experiments. 



ported data and can even alter the conformation of the 
macromolecule of interest. A theoretical understanding 
of these effects will improve the ability to use the TPM 
for quantitative |l3j analysis of biomolecular properties 
at the single molecule level. 

In the remainder of the paper, we first describe a simple 
statistical-mechanical theory of bead-induced volume- 
exclusion forces. We show how these forces depend both 
upon bead size and on tether length. We also derive 
scaling relations between the experimental mcasurablcs 
(bead position) and parameters such as bead size and 
tether length. Because the simple analytic model ne- 
glects some features of the full problem, we then turn to 
simulation results which capture all of the key effects and 
compare to the analytic results. 

The aim of the calculations outlined below is to il- 
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lustrate how the presence of a bead alters the statis- 
tical properties of the molecule to which it is tethered 
and how the bead reports information to the experi- 
menter. We confine our discussion to the equilibrium 
characteristics of this phenomenon, a key measurable in 
TPM experiments even for the study of dynamical pro- 
cesses |3 II II II II H EH We first note that in many 
experiments, the bead is flexibly linked to the end of 
its molecular tether, and hence is nearly free to rotate 
around the point of linkage 0, 0|. However, steric con- 
straints limit this freedom. In particular, the closer the 
bead is to the surface, the fewer angular conformations 
are available to it (Fig. [JJ . 

The statistical properties of the bead in TPM are de- 
termined by the coarse-grained free energy function (or 
"Hamiltonian" ) 

H = H m ({X}) + H b!Tn (R, {X}) + H b (K, {X}). (1) 

Here {X} is an abstract set of coordinates describing 
the configuration of the molecule and R is the vector 
pointing from the end-point of the molecule to the cen- 
ter of the bead (Fig. QJ. Eq.^contains three terms: The 
first describes the self-interactions of the molecule and in- 
teractions with external forces (surface forces or applied 
fields) other than those with the bead itself. Those in- 
teractions are captured by the second term, Hb. m - The 
last term, Hb, describes the external forces on the bead, 
for example those arising from applied fields or the sur- 
face. This term also depends on the configuration of the 
molecule, as the position of the bead depends on both its 
orientation R and the molecule's end-point. 

We obtain the statistical average of an observable A of 
the system as a weighted average over all the configura- 
tions of the system 
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Here / d 2 R is an integral over bead orientations and Z 
is the partition function. 

Before discussing the consequences of this model, we 
first discuss the significance of the terms in Eq.^ and 
make some initial simplifications. The external forces 
acting on the bead that we have in mind result from 
its interaction with the surface. This interaction con- 
tains the repulsive double layer pot ential and an attrac- 
tive van der Waals interaction [1J|, along with a hard- 
wall repulsion (Eq.El below) . Under physiological condi- 
tions, the double layer potential has an interaction range 
with a typical length scale of a nanometer, much shorter 
than the molecule lengths of interest to us; the van der 
Waals attraction, too, is weak on long scales |l5|. Ac- 
cordingly, we will model all bead-wall interactions using 
a simple hard-wall potential. We will also temporarily ig- 
nore bead-molecule interactions; later we will show that 
they have little influence for tethers like the ones of in- 
terest here. Given these assumptions, the last two terms 



in Eq.^ simplify to 
Hh m = 0, Hb 



if R(l - cos(a)) < z 
oo if i?(l — cos(a)) > z, 



(2) 



where z is the height of the end-point of the molecule 
and a is the polar angle of R (Fig. [TJ. 

We now examine the statistical averages of a molecular 
property A m ({X}) such as end-to-end distance. To ob- 
tain the average value of A m ({X}), not only do we need 
to sum over all of the configurations of the molecule, but 
also, we must sum over all of the configurations of the 
bead. Thanks to the simplifications in Eq.[3 the inte- 
gration over R can be done explicitly and results in an 
effective free energy function for the molecule. That is, 
the resulting statistical average of A m can be written as 

(A m ) = | J d{X}A m ({X})e-^«, where (3) 
fleff = H m ({X})-k B T\og(n(z)). (4) 

The new second term in H a f[ accounts for the configura- 
tions available to the bead. Q(z) is the solid angle allowed 
for R, given a molecular configuration {X}: 



Q(z) 



2-kz/R, z < 2R 
4tt, z > 2R. 



(5) 



The partition function is also consistent with the defini- 
tion of the effective Hamiltonian, Z = f d{X}e~@ eH . 

As a result of the constraints on the excursions of the 
bead there is an effective repulsive force, which prevents 
the end-point of the molecule from making contact with 
the surface and stretches the molecule. That is, the prob- 
lem is equivalent to one without the bead, but in which 
the end-point of the molecule is subjected to a force 
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where is the Heaviside step function. This entropic 
force alters the statistical properties of the tethered 
molecule (Eq.[3J, and can affect its interactions with itself 
or with other molecules. 

The key measurable associated with current TPM ex- 
periments is the position of the bead itself. That is, the 
output of the experiment is a record of the positions of 
the bead on successive video frames 0, UM LL2t UM ■ Hav- 
ing revealed that the confinement of the bead subjects 
the molecule to an entropic force, we now determine how 
this confinement influences bead motion. Note that one 
of our key conclusions is that there is a subtle dependence 
of the measured bead excursions on the size of the bead 
itself. To see this, let r (Fig. ^) denote the vector from 
the wall attachment point to the bead center. Given our 
simplifications (Eq.|2J), the non-zero statistical moments 
of r, up to second order, are 



(r z ) = (z)-±(z) e +R(l) e (7a) 
(r 2 ) = (^)_|^2 )0+jR ^ )e + i2 2 ( i + | (1)e) (?b) 

R(z) e 
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where rj_ is the in-plane displacement of the bead center 
and (xj_,z) is the displacement of the end-point of the 
molecule. Here all quantities averaged on the right-hand 
side are independent of the coordinates of the bead and, 
hence, are defined as in Eq.|31 The averages with sub- 
script correspond to summing over states only when 
the end-point satisfies z < 2R, i.e. 



(A 



m/e 



d{X}A m ({X})e-^Q(2R - z). (8) 



These 0-weighted terms arise because of the reduction 
of the configurations available to the bead due to the 
proximity of the surface (Eq.[5jl- Below it is shown how 
they give rise to experimentally testable scaling relations 
relating the excursion of the bead to its radius and the 
contour length of the molecule. 

We further pursue an analytic description of TPM by 
modeling the molecule as a Gaussian chain. The Gaus- 
sian chain is a useful approximation for molecules with 
short persistence lengths (such as RNA) and we demon- 
strate, using numerical simulations, that it also serves as 
a good guide for semi-flexible molecules {e.g. DNA) in 
the regime of interest here. 

Our problem involves a molecule grafted onto a surface, 
i.e. confined to a half space. DiMarzio showed that the 
corresponding Gaussian chain is described by 



(3fc B T/4L0(xx 2 + z 2 ) - k B Tlog(z/£). (9) 



L is the contour length of the molecule, £ is the persis- 
tence length and I is an arbitrary constant. The mate- 
rial properties of the molecule enter H m only through the 
combination L£. 

Because there are only two relevant length scales in the 
problem (R and y/L^) , we may write each moment of the 
molecule's excursions in terms of a function of a single 
variable. Using Eq.[!|]to evaluate the averages in Eqs.0 
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The excursions depend on the dimcnsionlcss number 
N-r, = R/ which we call the "excursion num- 
ber." Nn controls the bead's scaling behavior, defin- 
ing regimes of molecule-dominated motion (AT R < 1) and 
bead-dominated motion (JVr > 1, confined rotations). 

Figure [2] shows the relationship between excursions 
and iV R . For small excursion number, the excursions 
scale as (r 2 2 ) ps (r^ 2 ) ps L£ and (r z ) ps y/L^] the 
dependence on contour length obeys the expected rela- 
tions for a Gaussian chain with no bead attached. The 
scaling changes for large excursion number. Now, the 
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R/(L^/3) 1/2 

FIG. 2: Scaling behavior of bead excursion, normalized 
by coil size parameter, versus the excursion number Nn- 
Curves: analytical theory in the Gaussian-chain approxima- 
tion (Eg . 1101 . Circles: Monte Carlo calculation for a semiflex- 
ible chain with £ = 50 nm, L = 1245 bp, and various values of 
R. 
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FIG. 3: Solid curves: Theoretical prediction of the probabil- 
ity distributions for the projected distance r±, taking bead 
radius R = 250 nm, persistence length £ = 50 nm, and con- 
tour length L = 1000 bp (left curve) and 2000 (right curve). 
Dashed curve: Two dimensional Gaussian distribution with 
the same mean-square excursion as the left curve. 



mean excursions display power-law dependences on 7V R : 
(r ± 2 ) ps RVI4, (r z 2 ) ps R 2 and (r z ) ps R; the observed 
motion is dominated by the bead's rotation. The power 
law difference on bead radius R follows directly from 
the general relations (Eq.[?I), independent of the model 
representing the molecule. Now, the in-plane excursions 
show a square-root dependence on the contour length, in 
contrast with the linear dependence for small excursion 
number. This functional form arises because the average 
height of the molecule (which depends on y/L£) dictates 
the degree to which the bead can rotate. These scal- 
ing relations should be testable in experimental studies, 
where excursion numbers have ranged over wide set of 
values, 0.1 < N R < 90 0, II II H S IlS E3 (Most are in 
the regime 7V R > 1.) 

Confinement effects can alter molecular proper- 
ties, causing out-of-plane stretching of the molecule. 
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Eqs.H and H yield: 



L£/3 
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4 N R e~ N n 
7^ erf(JV R ) 



(11) 



In contrast, a Gaussian chain in free solution has 
squared out-of-plane excursions of the molecule equal to 
(£}/(L£/3) = 2. For many experiments H S H i, H ES 
lllj the excursion number is such that the second term 
in Eq.^1 is negligible and the effects of bead exclusion 
(Eq.|SJ) and molecule exclusion (second term in Eq.[3J) re- 
sult in a tripling of the out-of-plane squared displace- 
ment of the molecule. (Both exclusion effects contribute 
equally.) The bead- induced stretching can be viewed as 
a consequence of the effective force (Eq.[5| . In the Gaus- 
sian chain model, this force is 
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This force can significantly affect rates of loop formation 
in DNA. Finzi and Gelles used TPM to observe loop 
formation in DNA generated by the lac-repressor protein. 
Under their conditions, we predict an average effective 
force (F c g) 25 fN. Using the simple approximation that 
the rate of loop formation decreases by exp(— f3(F e g)l) (I 
is the operator-operator distance) we estimate that the 
bead-confinement effect reduces the rate of loop forma- 
tion by a factor of 2. 

To check the validity of our simplified model, we per- 
formed a simple Monte Carlo calculation. Our results 
agree with an independent calculation by D. Brogioli 
|l2j . Our code generated sets of discrete chains with ran- 
dom bends chosen to obtain a desired persistence length 
£. Each chain began at a random angle relative to the 
wall, and ended with the bead at a random orientation. 
Configurations where the bead, wall, or chain overlapped 
were discarded, (this includes bead-molecule interactions 
(Hb tTn ^ 0)) and the required averages were computed. 
Figure |21 shows that even for a stiff polymer like DNA, 
the scaling relations predicted by the approximate ana- 
lytical theory are accurate in the regime of interest to us. 
Actual experimental data allow the calculation of more 
subtle metrics than just averages, however: Fig. |31 shows 
predictions for the full probability distribution of excur- 
sions. The distribution is quite different from a Gaussian, 
a fact already observed experimentally [Tlj . 

The above results are interesting as fundamental poly- 
mer physics. For example, single-particle tracking al- 
lows the observation of the full probability distribution, 
and hence the opportunity to directly observe an end- 
end distribution for a semiflexible polymer and compare 
to our predictions. But our main goal was to develop 
a theoretical framework which can bolster the quan- 
titative capabilities of the TPM — a relatively noninva- 
sive, single molecule probe. We revealed that the prox- 
imity of the bead to the surface provokes an effective 



force on the molecule, altering its statistical properties 
and influencing biomolecular interactions. In addition, 
we determined how the excursions of the bead are in- 
fluenced by experimental parameters such as bead size 
and contour lengths; relations which are currently being 
tested 0. Finally, understanding how the competition 
between bead and tether effects is controlled by the ex- 
cursion number JVr, may help in the choice of optimal 
bead size and tether length for a given experiment. 
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